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We investigate the possible form of ideal intersections for two-dimensional rf trap networks suitable 
for quantum information processing with trapped ions. We show that the lowest order multipole 
component of the rf field that can contribute to an ideal intersection is a hexapole term uniquely 
determined by the tangents of the intersecting paths. The corresponding ponderomotive potential 
does not provide any confinement perpendicular to the paths if these intersect at right angles, 
indicating that ideal right-angle X-intersections are impossible to achieve with hexapole fields. Based 
on this result, we propose an implementation of an ideal oblique-X intersection using a three- 
dimensional electrode structure. 
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PACS numbers: 37.10.Gh, 37.10.Ty,41.20.Cv 

Intersections between network paths are a key ingredi- 
ent in two-dimensional (2D) rf trap networks which have 
been proposed to allow large scale quantum information 
processing (QIP) with trapped ions [1, 2, 3]. RF traps 
confine ions by a combination of rf and quasistatic electric 
fields [4, 5] and in an ideal network the rf field and the as- 
sociated ponderomotive potential vanish in the trapping 
zones and on dedicated paths between zones but nowhere 
else. Such ideal trap networks have been demonstrated in 
one dimension (ID) with segmented linear rf traps where 
transport as well as splitting and joining of groups of ions 
have been demonstrated to be possible with a very high 
degree of control [6, 7, 8]. In contrast, no ideal 2D trap 
networks have been identified, although a number of pos- 
sible intersection geometries for 2D trap networks have 
been investigated [9, 10, 11, 12, 13, 14]. The observed 
shortcomings fall in two broad categories. For T [9], Y 
[13, 15] and some X [10, 11], intersections, a residual 
rf field is observed in the paths through the intersection 
which can possibly lead to motional heating [16], in addi- 
tion to complicating or hindering controlled ion transport 
[12]. In contrast, some high-symmetry X geometries of- 
fer truly field-free paths but fail to be "fully confining" 
in the sense that there are unwanted lines of zero field 
allowing ions to escape the trap network [11]. 

In this paper we show that the unique simplest form of 
an ideal intersection for 2D trap networks is an oblique 
X and propose an implementation of this intersection 
based on a three-dimensional (3D) electrode structure 
that faithfully implements the ideal intersection. 

The paper is structured as follows: Sec. I defines the 
basic problem of designing ideal rf trap intersections. 
In Sec. II we constructively identify a hexapole term 
uniquely determined by the intersection angle as the only 
multipole term of hexapole or lower order which can con- 
tribute to a zero-field intersection. In Sec. Ill we investi- 
gate the properties of the identified hexapole intersection 
and show that fully confining hexapole intersections are 
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only possible at oblique intersection angles. Lastly, in 
Sec. IV we describe an implementation of the ideal inter- 
section based on a 3D electrode structure. 



I. RF TRAP NETWORKS 

Although rf traps use a combination of quasistatic and 
rf electrical fields to confine charged particles, we will 
here only be concerned with the rf component. For typ- 
ical QIP ion traps, the dimensions (< mm) and rf fre- 
quencies (> MHz) are such that we can adequately treat 
the rf field as quasistatic and express the time dependent 
rf field as cos(f2t) E(r), where f2 and E denote the rf 
frequency and spatially varying amplitude, respectively. 
The rf period 27r/f2 is made to be the fastest time scale 
of the trap, so that in the adiabatic limit the effect of 
the field can be described by an external ponderomotive 
potential U given by 
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where M and Q denote the mass and charge of the 
trapped particle [17, 18]. For the purpose of this pa- 
per, it suffices to note that U(r) is proportional to the 
square of the local rf field amplitude. Since confinement 
by electrostatic fields is impossible, the rf field must al- 
ways be present to allow trapping, and we will assume 
E(r) to be constant as is usually the case in QIP ion 
traps. In this case, U constitutes a constant "landscape" 
in which ions can be shuffled around by manipulating a 
quasistatic control field provided by dedicated trap elec- 
trodes. As mentioned above, we are here only interested 
in the structure of this landscape, which is determined 
entirely by the electrode geometry and is independent of 
the amplitude and frequency of the applied rf field. 

For an ideal trap network, as introduced in the In- 
troduction, we require that the E(r) vanishes for r on 
a network path and only there. Ideal networks have a 
number of qualities. Firstly, the ions do not experience 
any rf-induced micromotion. Micromotion is especially 
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critical in trap regions where gate operations are per- 
formed [19], but could also potentially lead to heating 
effects [16]. Secondly, the ponderomotive potential asso- 
ciated with a nonzero rf field along trap paths compli- 
cates fully controlled ion transport because the applied 
control fields must be engineered to compensate for any 
curvature of U along the path in order to avoid motional 
heating of the ions [8, 12]. Note that such compensation 
is impossible for the transport of multispecies crystals, 
which according to Eq. (1) experience different pondero- 
motive potentials due to mass differences. Lastly, the re- 
quirement that U(r) vanishes only on the network paths 
ensures that ions are confined, even in the absence of 
control fields. 

An ideal ID trap network can be implemented with 
segmented linear rf traps [6] where the rf field forms a 
quadrupole with a line of zero rf field along the trap axis 
[20], so that the isosurfaces of U are concentric cylinders 
around the axis. In this case, transport as well as split- 
ting and joining of groups of ions have been demonstrated 
to be possible with a very high degree of control [6, 7, 8]. 

We should note that although our quest in the fol- 
lowing is for ideal trap networks, none of the conditions 
for a network to be ideal in the sense introduced above 
are absolute requirements for a network to be suitable 
for even large scale QIP. Firstly, it has been proposed 
that transport in a network that is not fully confining 
can be achieved by "surfing" around any leaks [12]. Sec- 
ondly even in the presence of residual a rf field along 
the network path, it can be possible to perform highly 
controlled ion transport so that no or minimal heating 
takes place [8, 12]. Nevertheless, it seems that an ideal 
network would have a number of operational benefits. 

II. UNIQUE HEXAPOLE INTERSECTION 

As a first step to constructing an ideal 2D intersection, 
we will study the possible local structure of E(r) at the 
intersection of two smooth network paths in terms of a 
multipole expansion. Our goal is to find the lowest order 
multipole term that can contribute to E(r). Identifying 
the lowest possible multipole order is desirable for two 
reasons. Firstly, the confinement that can be produced 
for a given ion-electrode distance is stronger for lower- 
order multipole fields [15]. Keeping the ion-electrode dis- 
tance as large as possible is important since ion heating 
has been repeatedly demonstrated to increase dramati- 
cally with decreasing ion-electrode distance [21, 22, 23]. 
Secondly, the construction of an intersection where all 
low-order multipole terms vanish is technically challeng- 
ing. By trapping with the lowest order multipole possi- 
ble, a smaller number of terms must be made to vanish. 

With this in mind, we expand E(r) around the inter- 
section point chosen as the origin, 

E t = di + q itj rj + !/>,.,.;,/% n. + 0(r 3 ), (2) 

where the tensors q^j = djEi and hij.k — djdkEi de- 



scribe the quadrupole and hexapole components of the 
field. Here, as in the remainder of the paper, we have 
adopted the convention that summation over repeated 
indices is implied, e.g. qtjrj = q it ix + q^y + q i:3 z. 

It follows from Maxwell's equations that E(r) must 
be irrotational and divergence-free in the trap region, i.e., 
V x E(r) = and V-E(r) = 0. These requirements give 
rise to constraints on the quadrupole and hexapole ten- 
sors. To find these constraints, we note that the require- 
ments are equivalent to requiring that E(r) = -Vy(r) 
for an electric potential V fulfilling the Laplace condition 
didiV(r) = 0. For the quadrupole terms, this implies 
that qij must be symmetric and traceless, 

9i,j = Qj,i and = 0. (3) 

This leaves five free parameters as expected since a ba- 
sis of the quadrupole field components can be formed 
from the five second-order spherical harmonics. For the 
hexapole terms, we find similarly that hij.k must be sym- 
metric, and that the partial traces must vanish since 
didjdjVir) = d^Vir), 

hi,j,k = hi t kj = hj^^k and hijj = 0. (4) 

Here, the symmetry requirement leaves ten independent 
parameters with i < j < k, so that after the three partial 
trace constraints are included, seven free parameters are 
left as expected. 

We will describe the desired paths of zero field 
through the intersection by two curves 7^(s), I = 1,2, 
parametrized by path length s. The requirements that 
the curves intersect at the origin and that E vanishes on 
the paths can then be expressed as 

7 (0 (0)=O, (5a) 
E(~fW(s)) = 0,for all s and I = 1,2. (5b) 

To establish the resulting constraints on the multipole 
coefficients of E(r), we will expand 2£(7^(s)) in s and 
look at each order in turn. We will assume that the 
curves have a parametrization with continuous second- 
order derivatives. In this case, we can require the curves 
to be parametrized by path length, so that |9 s 7^'(s)| = 
1. Also, according to Eq. (5a) we then have that 

7 (0( s ) = f O s + N^k (1 Ks 2 + G(s 3 ), (6) 

where, and are the tangent and normal vectors 
of at the origin, and is the curvature at the same 
point. Except for the requirement that the tangent and 
normal vectors are perpendicular unit vectors, T^ l \ 
and «W can be chosen freely. 

Let us now assume that E(-y( L \s)) = as required by 
Eq. (5b) and deduce the possible value of the multipole 
expansion coefficients appearing in Eq. (2). Inserting 7^ 
in the form Eq. (6) into the expansion of E as given by 
Eq. (2), we find that 

E i { 1 «\s)) = d i + 0{s). (7) 
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For E( / y^ l \s)) to vanish to zeroth order in s, we must 
consequently have that 

(k = 0, (8) 

so that the homogeneous field component at the origin 
vanishes. Continuing the expansion of Ei(~f( l \s)), given 
that the d t fulfill Eq. (8), we find 

E i ( 7 U(s))=q i jfPs + 0(s 2 ). (9) 

It follows that for 2£( 7 W(s)) to vanish to first order in s, 
we must require that is an eigenvector of q im j with 
eigenvalue for I = 1,2. In the case where and 
T^ 2 ) are linearly independent (the paths are not cotan- 
gential at the intersection) , it follows from the additional 
constraint that qtj is symmetric and traceless [Eq. (3)] 
that 

Qi,j = 0. (10) 

A more detailed argument, presented in the Appendix, 
shows that the quadrupole components must vanish even 
in the cotangential case provided the paths are not iden- 
tical. For now, we assume the dipole and quadrupole 
terms to vanish as required by Eqs. (8) and (10), so that 
the expansion of £,(7' , '(s)) to second order takes the 
form 

^( 7 (O 00) - §/^, fe Tf f fs 2 + 0{s 3 ). (11) 

Note that the expansion only depends on the tangent 
vectors of the curves. For algebraic simplicity, we 
will orient the coordinate system so that 

f (1) = cos(0) x + sin(0) y (12a) 
f (2) = cos(0) x - sin(0) y, (12b) 

where 9 is half the angle between the tangents of the 
paths at the origin. Without loss of generality we require 
< 9 < 7r/4. In this case, we find by Eq. (11) that 

^( 7 «( S ))-^( 7 ( 2 )( S )) = 

2h iX2 cos(0) sin(0)s 2 + 0(s 3 ), (13) 

and consequently require = for £ l j( 7 ^) to vanish 
to second order. Similarly, we find that 

E i (^\ S ))+E i { 1 ^( S )) = 

cos 2 (0) + h it2 ,2 sin 2 (0)] s 2 + 0(s 3 ), (14) 

so that since by Eqs. (4) and (13) fti.,1,2 = ^1,2.2 = also 
/i2,2,2 and must vanish. Taking the partial trace 

constraint Eq. (4) into account, we find that the only 
hexapole terms that can be nonzero are /ii 1 3, ^-2,2,3, an d 
^3,3,3- For these terms we have, according to Eqs. (4) and 



(14), that E(j( l \s)) vanishes to second order in s if and 
only if 

c °s 2 W + ^2,2,3 sin 2 (0) = (15a) 
^1,1,3 + ^2,2,3 + ^3,3,3 = 0, (15b) 

with the unique solution that 

hi,i, 3 \ ( -sin 2 (0) \ 

^2,2.3 =6 a x +cos 2 (0) , (16) 
^3,3,3/ \sin 2 (#)-cos 2 (0)/ 

for some constant ax- To fulfill the requirements of 
Eqs. (8), (10), and (16), we find that the electrical po- 
tential must be of the form 

y(r) = a x e x (r)+0(r 4 ), (17) 
where the basis function 6x is given by 

6 x (r) = sin 2 (0) z(3x 2 - z 2 ) - cos 2 (9) z(3y 2 - z 2 ) (18) 

= sin 2 (0) e ( *\r) - cos 2 (0) 6^ } (r), (19) 

where we have introduced the basis functions 6^ = 
z(3r 2 — z 2 ) for later reference. 

To summarize, we have constructively established that 
to form a zero-field intersection, E(r) must have van- 
ishing dipole and quadrupole terms. Furthermore, if the 
intersection is not cotangential, the only possible field 
configuration is E(r) = -«xV6x(f) + 0(r 3 ). 

The key assumption in the arguments presented above 
is that the curves 7 *^ have continuous second derivatives, 
so that the expansion Eq. (6) is valid. We assume it to 
be the case that all isolated lines of zero field will have 
this property, but even though the structure of zero-field 
points in free space electrostatic potentials, or equiva- 
lently critical points of harmonic functions, has been in- 
tensely studied, most results appear to pertain to 2D 
systems and we have not found any direct argument to 
support our assumption, which must consequently stand 
as a conjecture. On the other hand, we can say for cer- 
tain that if the field vanishes on an interval of a curve 
with analytic parametrization, e.g., a straight line seg- 
ment, then it must vanish everywhere on the analytic 
continuation of that line segment. This follows from the 
fact that the maps s — > E^j^ (s)) will in this case be 
analytic as long as 7 ^(s) is in free space. This clearly 
rules out the possibility of ideal straight-line Y or T in- 
tersections where a straight line segment with zero field 
would have to terminate at the intersection point. 

Finally, we should point out that a fully confining, 
zero-field, right-angle intersection does exist. According 
to the results above, for this case V(r) must be C(r 4 ), so 
that the resulting intersection has a number of disadvan- 
tages compared to hexapole intersections as discussed in 
the beginning of this section. By brute force search in the 
possible octupole terms we find that V(r) = a Qo( r ) f° r 

o (r) = z i -3(x 2 + y 2 )z 2 + 3x 2 y 2 (20) 




FIG. 1: (Color online) Isosurface of |V6 X | 2 for 9 X of 
Eq. (17) for 8 = tt/6 (a) and 7r/4 (b). The isosurface is 
truncated to a cube: according to Eq. (22) the zero-field lines 
extend to infinity. Note that there is no confinement along 
the z axis for 9 = 7r/4. Since all terms of 0x are third order 
in the coordinates, the isosurfaces of |V0 X | 2 are related by 
scaling about the origin. 



is an example of an octupole potential giving rise to an in- 
tersection with zero field on the x and y axes but nowhere 
else. 



III. PROPERTIES OF THE HEXAPOLE 
INTERSECTION 

In this section we investigate the properties of the 
hexapole intersection described by Eq. (17) in more de- 
tail, first with respect to confinement properties and then 
with respect to implementation considerations. 



A. Confinement properties 

To study the confinement properties of U cx |VOx| 2 
corresponding to the hexapole field described by 0x, we 
look for points where the corresponding field 

/ 2xz sin 2 (9) 

V6 x (r) = 3 -2yz cos 2 (0) 

\(z 2 - y 2 ) cos 2 ((9) - (z 2 - x 2 ) sin 2 (6l), 

(21) 

vanishes. In the general case, this is seen to happen along 
the two lines 



y = ±tan(0)x for z = 0, 



(22) 



which are exactly the tangent lines specified in Eq. (12). 
However, for the special case of 9 — 7r/4, VGx also van- 
ishes everywhere on the z axis, as illustrated in Fig. 1. 

To understand the transition to vanishing confinement 
in the z direction, we consider the strength of the (quar- 
tic) confinement provided by U along the y and z axes 
as follows: 



ive> 
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y 4 cos 4 (9) 
z 4 cos 2 (26») 



for x = z = 
for x = y = 0. 



(23) 



Although confinement strength decreases monotonously 
with increasing 9 for both axes, only the z confinement 
vanishes completely. 





FIG. 2: Field lines of a hexapole and a homogeneous field, 
together with the field lines of the superposition field which 
has two quadrupole poles. Reversing the sign of the homo- 
geneous field would shift the two quadrupoles to above and 
below the hexapole center. 



B. Implementation considerations 

One reason that low-order intersections are interesting 
is that in practical implementations, it is difficult to pas- 
sively ensure that all lower-order terms vanish exactly. 
For typical traps with nonadjustable electrode geometry 
and a single rf feed only the strength and frequency of the 
rf field can be changed once the trap is in its operational 
configuration. The geometry of the rf field, including the 
relative strength of the multipole terms at the intersec- 
tion point, is not adjustable for such traps. Since lower- 
order terms dominate at the intersection point, this could 
be critical to the performance of the intersection. We will 
consequently investigate which lower-order terms are al- 
lowed by symmetry and try to describe the perturbing 
effects of these on the hexapole intersection, Eq. (17). 

Let us first note that the hexapole intersection poten- 
tial described by Qx(f) is symmetric in x and y and 
antisymmetric in z, 

@x(z, V, z) = 6x(-a;, y, z) = 6 x (a;, -y, z), (24a) 
6 x (a;, y, z) = -Q x (x, y, -z). (24b) 

An analysis similar to that of Sec. II shows that these 
symmetry requirements alone constrain V to be of the 
form 



V(r) =a H z + a x x (r) + C(r 4 ) 



(25) 



for some constant an. Under the generous assumption 
that an implementation exactly obeys the symmetries, 
Eq. (24), the only possible deviation from the ideal inter- 
section given by Eq. (17) is consequently a homogeneous 
field in the z direction. To empirically design a zero-field 
intersection under this symmetry, it is thus sufficient to 
ensure that the field vanishes at the intersection point, 
so that «h = 0. 

In realistic implementations, any field component not 
explicitly disallowed by symmetry would most likely be 
nonzero. It is consequently of interest to investigate the 
properties of U corresponding to the form of V given 
by Eq. (25) for small nonzero values of qh- An intuitive 
understanding of the effects of a perturbing homogeneous 
field can be had by considering the form of VOx in the 
y-z plane. Here by Eq. (21), 



V9 x (r) = 3 





-2yz cos 2 (9) 
z 2 cos(26») - y 2 cos 2 (6>), 



for x = 0, (26) 
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a) an/ ax < 

y 



b) tt H /ax = 

y 




located at the four points (all sign combinations) 



c) an/ ax > 

y 



FIG. 3: (Color online) Intersection structure modifications in 
the presence of a homogeneous field component, as described 
by oh in Eq. (25). Left-hand column: Contour curves of 
|V(aH-z + Qx©x)| 2 on x-y and x-z cut planes. Circles mark 
the barriers with a position given by Eq. (29) for «h/qx < 0, 
triangles mark the field zeros i n the y-z p lane as described by 
Eq. (27). The length scale is ^|a H /3a x | and 9 = 7r/6. Right- 
hand column: Full 3D form of the next-to-outermost contour. 
Note that in all three cases the outermost contour is connected 
and confined transport between any of the intersection entries 
is possible although the height of the barriers (circles) differs 
as discussed in the text. 



describing a deformed in-plane hexapole field. As illus- 
trated by Fig. 2, a homogeneous field will split a hexapole 
field into two quadrupoles. For a field of the form given 
by Eq. (25), the positions of these quadrupole field zeros 
are 




±860(6*) y 
±y/sec(26)z 



for an /ax > 
for an /ax < 0. 



(27) 



The positions are seen to depend critically on the sign of 
an- For an/ ax > 0, the zeros are located on the j/-axis, 
while for an /ax < they are located on the z-axis. 

The global shape of the intersection described by 
V(r) = anz + ax0x(f ) is very different in these two 
cases as illustrated by Fig. 3. For an/ax < 0, we get 
a double-junction structure with two junctions on the x 
axis, connected by paths in the x-z plane. The field is 
zero on the paths in the x-y plane, but this is not the case 
for the paths joining the junctions where we find barriers 
with a maximum height of 




±Vcsc 2 (6») -sec 2 (6»)x ± sec ((9) z 



(29) 



This structure is similar to that described in Refs. [24, 
25]. For an/ax > 0, we get two disjoint paths through 
the intersection. The height of the barrier at the origin 



|W| 2 



(30) 



corresponding to the pondcromotive potential due to the 
homogeneous field component. 

Given a certain production tolerance, Eqs. (28) and 
(30) show that it is beneficial to aim for a negative value 
of an /ax, corresponding to the double-junction config- 
uration. For a given value of |o:h|, this will reduce the 
barrier height by a factor of tan 2 (ff) compared to the case 
of an /ax > 0. 



IV. IMPLEMENTATION WITH 3D 
ELECTRODE STRUCTURE 

We will now describe a 3D electrode structure that can 
be used to implement the ideal intersection described by 
Eq. (17). 

Firstly, we note that the full symmetry of Ox as de- 
scribed by Eq. (24) can be retained in an implementation 
based on a 3D configuration of electrodes as illustrated 
in Fig. 4. According to the results of the previous sec- 
tion, it follows that the potential at the origin must be 
of the form given by Eq. (25). In addition, it turns out, 
the symmetry ensures unbroken lines of zero field in the 
x-y plane. By Eq. (24) the field in the x-y plane will be 
along the z direction since the x and y components of 
the field will be odd in z. Furthermore, the z component 
of the field in the x-y plane is a continuous function of x 
and y and the zero contour of this function, consisting of 
unbroken lines, will define the points of zero field. 

In total, these properties ensure that it is relatively 
easy to design a zero-field intersection using 3D elec- 
trodes: After choosing an overall design that ensures that 
there are paths of zero field leaving the intersection point, 
all that remains is to ensure that an = —E z (0) = 0. 
In the implementation illustrated in Fig. 4, this was 
achieved by a bisecting search for the correct value of 
a single dimensional parameter. 

It should be noted that symmetry alone does not define 
the angle 8 between the intersecting paths, but rather 
constrains the hexapole component to be a superposition 
of two components 0^ (r) and 0^\r) as follows: 



V(r) = a ^\r) + a^\r). 



(31) 



|W| 2 = an tan 2 (6>), 



(28) 



For the implementation illustrated in Fig. 4, we numeri- 
cally find that a^ ] = 0.011 V Q /d 3 and a^ = 0.12 V /d 3 , 
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FIG. 4: (Color online) Example implementation of an ideal 
intersection using a 3D electrode structure, (a) 3D view of 
the electrodes (red and gray polygons) and isosurfaces of 
| VV(r)| 2 (blue to green curved surfaces) with all components 
with z > cut away, (b) Cuts through the electrode struc- 
ture, which is symmetric under inversion in the x-z and y-z 
planes, and antisymmetric under inversion in the x-y to main- 
tain the symmetry of 0x as described by Eq. (24). Units are 
the ion-electrode distance d. Results are based on numerical 
simulations performed using the CPO program [26]. 



where V$ is the peak amplitude of the applied rf and d is 
the distance from the trap center to the nearest electrode. 
These values correspond to the rf field being described 
by Eq. (17) with a a strength of ax = 0.13 Vo/d? and an 

angle of 6 = arctan(^/ <*x V a x^) = ■ The obtained 
strength seems reasonable, given that the strongest pos- 
sible hexapole guides along the x or y axes would cor- 
respond to a^"' = ^Vo/d 3 [15]. The angle is rela- 
tively close to the asymptotic angle of 9 = 15° cho- 
sen for the demonstration design. If better agreement 
was desired, an additional fit parameter (for example, 
the length of the bridge recession) could have been in- 
cluded to ensure that the unwanted hexapole component 
cos 2 (9)Q^ + sin 2 (9)0^ associated with a given value 
of 9 was made to vanish. 



V. CONCLUSIONS AND OUTLOOK 

In conclusion, we have shown that a zero-field inter- 
section for rf traps cannot have any field component of 



quadrupole or lower order at the intersection point. Fur- 
thermore, we have demonstrated that the hexapole field 
component at the intersection point is uniquely deter- 
mined by the intersection angle [Eq. (17)], and that this 
component does not provide any confinement perpendic- 
ular to the intersecting paths if these intersect at right 
angles. These results can serve as a guide for the de- 
sign of intersections for rf trap networks suitable for QIP 
based on trapped ions. We have suggested how the in- 
tersection could be implemented using 3D electrodes. In 
relation to implementations of the intersection, we have 
shown that if static field components cannot be guaran- 
teed to completely vanish at the intersection point, it is 
beneficial to ensure that the imperfect intersection is of 
the "double-junction" type [Fig. 3 (a)]. 

The proposed implementation of the ideal intersection 
illustrated in Fig. 4 is not well suited for microfabrication, 
which would most likely be a requirement for fabrication 
of large scale trap networks. For this reason it would 
be worthwhile to find intersection implementations with 
electrode geometries better suited for microfabrication, 
such as the surface-electrode geometry [11], which is com- 
patible with large scale microfabrication [27] , and has the 
lowest demonstrated heating rates for microfabricated 
traps [23, 28, 29]. Initial investigations in this direction 
indicate that such intersections need to be very oblique 
(9 as 7t/12) to obtain reasonable hexapole strengths for 
an intersection of two straight guides. As an alterna- 
tive, it might be worthwhile to investigate cotangential 
intersections of curved guides. 
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APPENDIX A: GENERAL PROOF OF ABSENCE 
OF QUADRUPOLE TERMS IN AN 
INTERSECTION 

Here we present a general proof that there can be no 
quadrupole terms in the multipole expansion of E as 
given by Eq. (2), even in the cotangential case where 
= T( 2 \ provided that the paths 7^ have analytic 
parametrizations and are not identical. 

We expand the paths as 

00 

7 W(s)=^wWs" (Ala) 

n=l 

00 

7< a >(«) =J2v (n) s n , (Alb) 

n=l 

where in particular — and t?W = T^ 2 ' are the 
identical tangent vectors. We now take M to be the 
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smallest integer so that tt( M ) ^ v^ M \ Such an M must 
exist since the paths are assumed not to be identical. 
Assume as in Sec. II that E(~f®(s)) = for all s so that 
by Eq. (8), d t = 0. We then find that 

q i , j (u'. M) -v^)s M + 0(s M +% (A2) 

so that for the difference to vanish to M-th order, we 
must have 



t (M) _ V (M) g N 



(A3) 



where N q is the null space of the q tensor. As in the 
noncotangential case, we have from Eq. (9) that 



t (D = V (D e N 



(A4) 



On the other hand, letting j(s) = Yl n <M u ^ sn we 



have that 



9s7 M( s ) =\d s j(s)f 



+ 2u« • tt^/- 1 + 0{s M ), (A5a) 
d sl ^(s)\ 2 = \d s j(s)\ 2 

+ 2vW ■ v^s M ~ l + 0(s M ). (A5b) 
Taking the difference, we find by Eq. (A4) 

d s7 (1) (s)| 2 -|d s7 (2) (s)| 2 = 

2«« • («( M > - s^- 1 + G(,s M ), (A6) 

so that for the difference to vanish, which it must since 
both curves are parametrized by path length, we must 
have 



"•( 



U {M) _ V (M) 



0, 



(A7) 



so that and u^ M ^— v( M ) are linearly independent. As 
and vS M > — «(*0 are both non-zero and by Eqs. (A3) 
and (A4) members of N q , it follows that qij = since 
qij is symmetric and traceless. 
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